Discontinuous conductance of bichromatically ac-gated quantum wires 
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We study the electron transport through a quantum wire under the influence of external time- 
dependent gate voltages. The wire is modelled by a tight-binding Hamiltonian for which we obtain 
the current from the corresponding transmission. The numerical evaluation of the dc current reveals 
that for bichromatic driving, the conductance depends sensitively on the commensurability of the 
driving frequencies. The current even possesses a discontinuous frequency dependence. Moreover, 
we find that the conductance as a function of the wire length oscillates with a period that depends 
on the ratio between the driving frequencies. 

PACS numbers: 05.60.Gg, 73.23.-b, 73.63.Nm 



I. INTRODUCTION 

Recent advancements in nanotechnology enabled the 
investigation of low-dimensional systems like quantum 
dots, quantum wires or few-atom systems. Thereby, 
techniques like mechanically controlled break junc- 
tions, scanning tunneling microscopy, or angle-resolved 
photo-emission electron spectroscopy provide informa- 
tion about the electronic structure or transport proper- 
ties of these systems. Particularly interesting are the 
transport properties of low-dimensional nanoscale con- 
ductors in the presence of time-dependent external fields 
for which a wealth of effects has been predicted like, e.g., 
photon-assisted tunneling,—!^ and electron pumping.— 
External fields can be also used to control both the dc 
current^ and the zero-frequency noise^-^ in mesoscopic 
conductors. 

In undrivcn one-dimensional quantum wires, con- 
ductance oscillations — the oscillatory dependence of 
the conductance as a function of the wire length — 
were predicted^i^ and experimentally confirmed^iS, They 
emerge when the energy levels of a quantum wire match 
the Fermi energy of the system and relate to charge neu- 
trality of the wire.JJi The period of these oscillations 
is usually two sites for undriven wires, but may be- 
come larger if the wire is exposed to a time-dependent 
field^^^ 

A second harmonic added to an ac field that drives 
a quantum wire may significantly alter transport proper- 
ties. The main reason of the differences is that the second 
harmonic may brake time-reversal symmetry and gen- 
eralized parity dynamically, leading to effects like non- 
adiabatic electron pumping ^li^iiiii^ii^ i.e., to a dc cur- 
rent even in the absence of any net bias. However, if 
the ac field acts as a gate voltage that shifts all wire 
levels uniformly, the pump current must vanish.^® Sim- 
ilar effects have been predicted for quantum Brown- 
ian motion in bichromatically driven infinitely extended 
periodic potentials.— i^i^^ Moreover, two incommensu- 
rable rectangular-shaped driving forces can cause a dc 
curren t^^i^^ and lead to interaction-induced frequency 



mixing. Also dynamical localization effects in quan- 
tum dots may be affected by bichromatic driving fields, 
in particular when the freuquencies of the latter are 
commensurable<2^i21 

A well established theoretical tool for treating period- 
ically driven quantum system is Floquet theory, which 
is based on the discrete time-translation invariance of 
the driving field. Higher harmonics of the fundamental 
driving frequency do not affect this invariance and, thus, 
can be included straightforwardly. As soon as the ratio 
between the fundamental frequency and its higher har- 
monic increases, a numerical treatment may require quite 
some effort. The situation becomes even worse, when 
the frequencies are no longer commensurable. Then it 
is rather desirable, even from a practical point of view, 
to find analytical expressions for the behavior of the 
driven quantum system, as for example in the recent 
works 20,21,22. For the description of electron trans- 
port through monochromatically driven systems, com- 
mon transport theories have been combined with Floquet 
theorySi^s^ 

In this paper we study the transport properties of a 
bichromatically driven quantum wire, focussing on the 
influence of commensurability of the driving frequencies. 
In doing so, we generalize in Sec. |TT] recent theoretical ap- 
proaches such that they allow one to treat two ac fields 
with arbitrary frequency ratio. Then we investigate in 
Sec. lIIII the influence of bichromatic driving on the trans- 
mission and the conductance oscillations. 



II. THEORETICAL DESCRIPTION 
A. Wire-lead model 

We consider the time-dependent wire-lead Hamiltonian 

H{t) = H^{re(t) -I- iJleads + -^coupl, whcrC 
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H^irc{t) = ^[l^n + f{t)]cXcn+^ Kc+C„+i -|- h.C. (1) 
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models in a tight-binding description the wire with sites 
|n), where neighboring sites are coupled via tunnel matrix 
elements Vn- The on-site energies e„ are modulated by 
an ac gate voltage f{t). The first and the last site of the 
wire are coupled via the tunnel Hamiltonian 



coupl 



(2) 



kL 



kR 



to leads at the left (i) and at the right (i?) end. 
Within a wide-band approximation, we assume that 
the spectral density of the wire-lead coupling, Ta{e) = 
2^Sfc l^fcP'^(£ ~ ^a) = Tq, is energy independent. The 
leads are modelled as ideal Fermi gases by the Hamilto- 



leads— ^ ^ka'^ta'^ka- 
a=L,R g 



(3) 



The operators c^;^, c„, and their Hermitian adjoints are 
the usual annihilation and creation operators for an elec- 
tron in lead a — L, R with wave vector k and for an 
electron at the nth wire site, respectively. We focus 
on experiments like those described in Refs. [ol lsills^ . 
in which electron-electron interactions can be neglected. 
For photon-assisted transport, this is, e.g., the case when 
the photon energy is much smaller than both the typical 
splitting between single-particle energies of the wire and 
the Coulomb interactioni^^ 



B. Time-dependent scattering theory 

Following the scattering approach of Ref. [l6|, we derive 
from the time-dependent wire-lead Hamiltonian H(t) the 
Heisenberg equations of motion for the creation opera- 
tors. The wire part of their solution can be expressed in 
terms of a retarded Green function G(t,t') which, owing 
to the time-dependence of the Hamiltonian, depends ex- 
plicitly on both times and fulfills the equation of motion 



^n— -Ho- f(t) + -r 



G{t,t')^5{t-t'). (4) 



The Hamiltonian Hq = ffwire — /(i)lwirG describes the 
wire in the absence of the time-dependent gating, while 
r = ri|l)(l| + F/jI A^) (A^ I is a self-energy stemming from 
the coupling of the wire to the leads. 

The solution of Eq. ([4]) can be traced back to the solu- 
tion of the time-independent problem for which the Green 
function g obeys the equation of motion 



g{t-t') = 6{t-t'). 



(5) 



Then it is straightforward to show that the Green func- 
tions G and g are related according to 



(6) 



This means that the influence of the driving field is sub- 
sumed in the phase 



dsfis) 



(7) 



Defining the current through the wire as the change of 
the electron number in the left lead, / ~ —ednL{t)/dt, 
we obtain for its expectation value the time-dependent 
Landauer-like expression 

Ht) - de[TRLit,e)hie) -TRLit,e)h{e)] - QLit) 

(8) 

with the time-dependent transmission 



TLRit,e)=TLTn.\Gi,Nit,e)\^ 
and Tj^R defined accordingly. The function 



(9) 



G{t,s) ^ j dTe'"G{t,t-T) (10) 
= e-*^(*) / dTe*^^+'^(*-^)g(r) (11) 



is the Fourier transformed of the Green function with 
respect to the time difference r = t — t'. Note that G still 
depends explicitly on the final time t. The displacement 
current qL{t) describes charges oscillating between the 
left lead and the wire and does not contribute to the dc 
current considered below. Thus it will not be considered 
henceforth. 

In order to evaluate the r- integration in Eq. (jlip . we 
have to specify the time-dependent gate voltage f{t). We 
assume that it consists of K harmonics, such that 



K 



fit) = ^ Afecos(wfet), 
fe=i 
K A 

Pi^) = V— sin(c^fci), 



fe=i 



(12) 
(13) 



where LOk and A/j are the frequency and the amplitude of 
the fcth harmonic of the ac field. Then the Green function 
PT|) reads 



G(t,e)= ^ J,nA^l/uJl)---Jm^{^K/i^K) 



mi,--- ,nij,^ 



exp (i ■mktJkt)g{e - ^ mkUJk) ■ 
k 



(14) 



k 



The Bessel functions Jm stem from a Fourier decompo- 
sition of the exponential in Eq. ([HI), while the evaluation 
of the r integral provided the Fourier transformed of the 
Green function, g{e). Note that we have ignored a time- 
dependent but nevertheless irrelevant phase factor. 

Assuming moreover symmetric wire- lead coupling, i.e., 
Fi = T n = F, the time-dependent transmission ^ be- 
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comes 



integers. Then we obtain 



mi ,■ ■ ■ ,mK 

k k 

X exp (i ^(mfc - ■m'^)ujkt) (15) 
= Tfli(t,e). (16) 

The relation in the last line means that from the right 
lead to the left lead has the same probability as the re- 
versed process. This is particular for ac gating and does 
not hold for general driving.— Here it stems from the 
fact that for the present model with time-dependent gat- 
ing, the time-dependent transmission can be expressed in 
terms of the Green function g for the time-independent 
problem. 

The dc current, being the main quantity of interest, is 
then obtained from Eq. ([5]) and reads 



de[fLie) ~ fLie)]T{e 



(17) 



where T{e) denotes the time-average of the transmis- 
sion Tifl(t,e), and its computation from Eq. (|15p seems 
straightforward. This is however not the case, because 
the averaging procedure depends crucially on the com- 
mensurability of the frequencies tok- If the gate volt- 
age f{t) consists of many spectral components, the ac- 
tual computation can be quite cumbersome. However, 
we below restrict ourselves to two frequencies, so that it 
is sufficient to address two cases: all frequencies being 
commensurable or all being incommensurable. 

For incommensurable frequencies, the phase factor in 
the last line of the time-dependent transmission (jl5p van- 



ishes if and only if 



for each mode k. Conse- 



quently, we obtain the average transmission 



r(e) = T'JL{^i/^i)---JIA^k/ook) 



m 1 , ■ ■ ■ ,771/^' 



k 



(18) 



The Green function of the undriven wire, g(e), depends 
on the specific wire model and will be determined below. 

If the frequencies of the external fields are commensu- 
rable, time-averaging unfortunately does not result in a 
less involved expression for the transmission. Therefore 
we provide the explicit result only for the case of two fre- 
quencies with ratio loi/lo2 = m/m' , where m and m! are 



T(e)= rV„,(Ai/wi)J„,(A2/w2) 

<mi ,7712 ,m2> 

X J„,.(Ai/u;i)J™- (Aa/t^a) (19) 
X 9in{£ - miUJi - 777,2^2) 
X 5i7v(£ - ™Wi - ^A^2), 

where the angular brackets restrict the sum to those val- 
ues of mi, 7712, 7773, TO4 that comply with the condition 
(777,1 — 777,'i)u;i = —(7772 — 7772)^2- Note that for commen- 
surable frequencies, this condition can also be fulfilled 
for 777fe ^ m'l^. Since for incommensurable frequencies 
the restriction to the sum is fulfilled only for 7771 ~ m'l 
and 7772 = '772, expression (|19p incidentally holds for both 
cases. 

The transmissions ^TS\\ and (fTO|) represent the main 
analytical results of this work. In the case of commensu- 
rable frequencies, the driving is periodic and, thus, can 
be treated within Floquet scattering theory.^'--'-- Nev- 
ertheless, the explicit evaluation of the transmission can 
be quite cumbersome. The case of incommensurable fre- 
quencies generally even requires a "two-color Floquet the- 
ory" which can be numerically demanding. Only the re- 
striction to models with uniform gating allowed us to 
obtain the present compact expressions for the transmis- 
sion. 



III. WIRE WITH EQUAL ON-SITE ENERGIES 

In order to present explicit results, we have to spec- 
ify the wire model. A A^-site model that allows even 
further analytical treatment is one with equal onsite en- 
ergies £„ = £0 for all sites |77), 77 = 1, . . . ,N and equal 
couplings Vn — V. For such a wire, the relevant element 
of the Green function in the absence of driving reads 



5iAf(£) 



(-1) 



Af-l 



Vun{x) + iTuN-i{x) - T'^U]y^2{x)/4:V' 



(20) 

where x = {e — eo)/2V, while m„ denotes the 77th 
Chebyschev polynomial of the second kind. For a deriva- 
tion, we refer the reader to Ref. [S^. This relation together 
with Eqs. (HH]) and allows one to obtain explicit ex- 
pressions for the transmission of the ac-gated wire p]). 

In our calculations we use the effective wire-lead cou- 
pling strength F as energy unit. Then for F = O.leV, 
the current unit becomes 2eF/fi ~ 50 /iA. Moreover, we 
measure level energies with respect to the Fermi energy, 
which implies Ep = 0. Henceforth, we present results for 
a bichromatic driving field. 
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FIG. 1: Transmission for a wire with N — 5 sites with equal 
onsite energies eo = as a function of the frequency uji for 
the driving amphtudes A2 — 4 (a) and A2 = 2 (b). The 
tunnel coupling between two neighboring wire sites is V = 4, 
while the other driving parameters read UJ2 = 2, Ai =4. 
Solid lines are obtained from expression II 8 1 which is valid for 
incommensurable frequencies only, while the dashed lines and 
the circles are obtained from Eq. (|19|l and, thus, are valid for 
all driving frequencies uji, uJ2- The fractions in panel (b) mark 
selected frequency ratios 1^)2/(^1. 



FIG. 2; Average current for the chemical potentials fiL = 
—fJ'R = 0.1 as a function of the equal onsite energies £0 for a 
wire with A*' = 3 sites and tunneling matrix elements V — 4. 
The first driving field is kept constant [Ai = 4 with frequency 
LJi = 2 (a) and tJi = 4 (b)]. The second field has driving am- 
plitude A2 — 4, while its frequency is 0^2 = 4 (thick solid lines) 
and U12 = 3.99 (thin sohd lines), respectively. The dashed line 
in panel (a) marks the undriven case which corresponds to 
vanishing driving amplitudes, i.e. Ai = A2 = 0. 



A. Discontinuous transmission and corresponding 
current 

In order to reveal the role of commensurability of the 
driving frequencies, we investigate the transmission of a 
wire with onsite energies at the Fermi energy level. Fig- 
ure[T]shows the result as a function of the frequency of the 
first external field, wi, while the frequency of the second 
field is kept at a constant value. As a most prominent 
feature, the transmission turns out to be a discontinu- 
ous function, which assumes exceptional values whenever 
LOijujx is rational, i.e., when the driving frequencies are 
commensurable. A clearly visible effect is found for ra- 
tios that can be expressed by small natural numbers. In 
particular, for the ratios 1/1 and 2/1, the deviation from 
the transmission in the vicinity of the resonance is of the 
order of 20%. For "less rational" numbers like e.g. 4/3 
or 2/3, we still find a small discontinuity which, however, 
is hardly visible on the scale chosen. The values can de- 
pend sensitively on details like driving amplitudes and 
wire parameters. In the present case, it is even such that 
the deviation at the prime resonance 1/1 turns with in- 
creasing amplitude A2 from a remarkable negative value 
to a positive value; cf. panels (a) and (b) of Fig. [T] More- 



over, the 2/1 resonance, albeit "very rational" , does here 
not lead to a significant discontinuity. A similar behavior 
has been observed for the transport in classical ratchet 
devices It is worth noting that the transmission for 
incommensurable frequencies [Eq. (|18p and solid lines in 
Fig.[T] possesses local maxima at the resonances 1/2, 1/1, 
and 2/1. 

In order to investigate further the role of commensura- 
bility, we focus on the 2/1 resonance for ijJ\ — 2 and the 
1/1 resonance for lo\ = 4. In doing so, we compare the 
current for the commensurable case 0^2 = 4 with the one 
found for the "less commensurable" uj^ — 3.99. The cur- 
rent for a wire of length = 3 as a function of the onsite 
energy eg is shown in Fig. [21 A common feature of both 
cases is the emergence of side peaks: The undriven wire 
possesses the eigenenergies Eq = Eq and E± = Eq ± 5.6. 
Whenever these eigenenergies lie within the voltage win- 
dow at the Fermi energy Ep = 0, i.e. for — Q and for 
£0 = T5.6, the current assumes a maximum in compli- 
ance with the Breit-Wigner formula; see Fig. 1 of Ref . [34l. 

In the presence of ac gating, we find side peaks shifted 
from the original peaks by multiples of the driving fre- 
quencies LUi and L02- Overlapping peaks may lead to de- 
structive interference. This is, e.g., the case for the data 
shown in Fig. [DJa) in the vicinity of eg = 0, where the 
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side peaks at eq T T ^2 for eq = ±5.6 compensate the 
fundamental peak at eg = 0; such that the total current 
is characterized by the local minimum. The proper loca- 
tion of the side peaks is most clearly visible for uji — 4 
and ^2 ~ 4 [Fig. [U^a)], i.e., when both frequencies are 
practically indistinguishable. Then we observe side peaks 
at Eq = 4, 8 (main peak at Eq = 0), at Eq = 1.6, 9.6 
(main peak at Eq = 5.6), and at Eq = 2.4 (main peak at 
Eo = —5.6). Let us emphasize that the location of the 
peaks does not allow us to divide the commensurable case 
from the incommensurable case, because the resonance 
width r = 1 is much larger than the difference between 
the two values for uj2 which reads 4 — 3.99 = 0.01. 

Despite the practically identical location of the peaks, 
the peak heights are significantly different for the two 
cases. This corresponds to the discontinuity observed al- 
ready in Fig. [1] The difference can even be up to a factor 
5, as is the case for wi = 4 [Fig.[2I^b)] for the zeroth-order 
resonances at eq = and Eo = 5.6. This means that in 
a possible experiment, the current may be controlled via 
a slight shift of one of the driving frequencies, provided 
that one tunes the wire parameters such that the current 
in the absence of the driving assumes a large value. 



B. Conductance oscillations 

The conductance G of a tight-binding system with 
onsite energies at the Fermi level obeys oscillations as 
a function of the wire length. Typically these os- 
cillations have period 2 (even-odd oscillations) 
but also larger periods have been predicted for both 
statio^i2£i2£i22 and monochromatically driven^^ wires. 
We here investigate the corresponding behavior for 
bichromatic gating. 

The oscillation period for the driven case can be esti- 
mated along the lines followed in Ref. [H: For a static 
wire, the emergence of conductance oscillation with pe- 
riod M means that a wires with N sites possesses the 
same conductance as a wire with N + M sites, i.e. 

5i,Af(e) = 5i,Af+M(e)- (21) 

After some algebra, one obtains^i cos(m7r/Af) — {e — 
Eo)/21/, where m = 0, 1, . . . , 7\f — 1. The energy differ- 
ence on the right-hand side of this equation marks the 
distance of the onsite energies Eq to the Fermi energy of 
the leads attached. In the presence of ac fields, onsite 
energies acquire shifts by multiples of the driving field 
quanta. Then the oscillation period is determined by the 
condition 

/m7r\ £ - + kiUJi + k2UJ2 

'=°nAF) = w ' 

where again m = 0, 1, . . . , Af — 1, while fci and k2 are 
integer numbers. 

In order to corroborate this condition, we have com- 
puted the transmission [which relates to the conductance 
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FIG. 3: Transmission as a function of the wire length A'' for 
Ai — A2 = 6, eo = 0, = 4. (a) the driving frequencies 
are loi = V^V, uj2 ~ 2V (circles) and ui — V , UJ2 = 2V 
(triangles) correspond to the oscillation periods Af = 4 and 
M = 6, respectively, (b) ui = 2V^ cos(7r/5), ijJ2 ~ '2V (circles) 
and uj\ — V , LU2 = \/2V (triangles) corresponding to Af = 10 
and M = 12, respectively. The lines serve as a guide to the 
eye. 

according to G = (2e^//i)r] for bichromatic driving. The 
results for various combinations of the driving frequen- 
cies are shown in Fig. [3l For all parameters used in the 
numerical evaluation of the transmission (|19p , we observe 
an oscillation period Af > 2 in compliance with Eq. (|22p . 
Conductance oscillations with period Af = 6 have also 
been found in Ref. [13 for wires under the influence of a 
monochromatic ac field. There the oscillation period is 
also determined by Eq. for /c2 = 0. We like to em- 
phasize that in the absence of the driving, by contrast, 
the oscillation period is Af = 2, i.e., we find the usual 
even-odd oscillations. 

A remaining intriguing question is whether for bichro- 
matic driving, the commcnsurability of the frequencies 
influences the conductance oscillations. Therefore, we 
compare for loi = V the commensurable uj2 — oJi 
(which in fact represents monochromatic driving) with 
the "practically incommensurable" UJ2 = 0.990)1. For 
both cases, Eq. (|^^ predicts conductance oscillations 
with a period of approximately Af = 6. Figure |4] de- 
picts the resulting transmission as a function of the wire 
length. Note, that the oscillation period for the static 
conductor is Af — 2, while it becomes Af = 6 in the 
presence of the driving. Generally, the variation of the 
conductance within one oscillation period is larger for 
incommensurate driving frequencies. For the driving pa- 
rameters used in Fig. 31 it is even such that for the com- 
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FIG. 4: Transmission as a function of the wire length A'' for 
commensurable frequencies {uii — ui2 = 4, triangles) and in- 
commensurable frequencies {uj2 ~ O.QQcui, uii = 4, circles) 
and the driving amplitudes Ai = A2 = A = 4 (a) and 
Ai = A2 = A = 8 (b). The other parameters are eo — 0, 
V — 4. The lines are a guide to the eye. 



mensurate lui =0^2, the oscillations have almost period 
M = 3 or M = 2 (solid lines) with a small variation of 
every second peak, such that strictly speaking, the pe- 
riod is still M = 6. Moreover, we again observe one 
significant difference between the commensurate and the 
incommensurate case: For given driving amplitudes, the 
conductance may change by almost one order of magni- 
tude upon shifting the second frequency from uj2 — ^^i to 
UJ2 = 0.99a;i. This again allows one in an experiment to 
switch the current by changing one driving frequency by 
only 1%. 

Note also that when the driving amplitude becomes 
smaller, the situation tends towards the undriven limit 
for which the oscillation period is M = 2. The transition 
is such that the conductance within the M = 6 period 
changes until eventually the oscillation period becomes 
practically M = 2. 



IV. CONCLUSIONS 

Using a Green's function technique, we have derived a 
Landauer-like formula for the current through a poly- 



chromatically ac-gated quantum wire described by a 
tight-binding model. If the gating acts uniformly on all 
wire sites, we have provided an expression for the time- 
dependent transmission that depends on the transmission 
in the absence of the driving and some time-dependent 
prefactors. The time-average of the transmission for ac 
gating determines the dc current being the central exper- 
imentally accessable quantity. In order to perform the 
time-average, we had to distinguish between commensu- 
rable and incommensurate frequencies. We have found 
that the transmission and, thus, the dc current depends 
sensitively on the commensurability of the driving fre- 
quencies. The difference is most noticeable close to the 
resonances with the "very rational" frequency ratios 1/2 
and 1/1, where the conductance as a function of one driv- 
ing frequency is discontinuous and may suddenly change 
by up to one order of magnitude. 

The commensurability of the driving frequencies leaves 
also its fingerprints in the conductance oscillations: 
While for static conductors, one usually observes an oscil- 
lation period of two sites (even-odd oscillations) , bichro- 
matic driving may increase the period significantly. For 
our model, we observed oscillation periods of up to 12 
sites. The physical reason for the enlarged period is that 
the bichromatic driving influences the resonance condi- 
tion that determines the oscillations. We have confirmed 
this intuitive prediction by numerical studies. As a re- 
markable influence of the commensurability, we found 
that the growth of the oscillation period is more visible 
for incommensurate frequencies. 

In a possible experiment with STM tips that contact a 
quantum wire on a vicinal surface ^^i^ the observed ef- 
fect may be used for switching the current upon slightly 
changing the frequency of the ac gate voltage. Alterna- 
tively, such an experiment may also be performed with 
a molecular wire that bridges a break junctio n^'^i'^^ or 
with a chain of coherently coupled quantum dotsi^i^ In 
particular the latter type of experiment is well suited for 
applying time-dependent gate voltages. 
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